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ABSTRACT:  Published  annual-to-decadal-resolution  stalagmite  81sO  records  since  AD  900  from  six  caves  (Dongge, 
Furong,  Heshang,  Buddha,  Shihua  and  Wanxiang)  in  China  were  analyzed  to  detect  temporal  and  spatial  variability  of 
the  East  Asian  Summer  Monsoon  strength,  which  strongly  affects  wet/dry  conditions  in  eastern  China.  The  empirical 
mode  decomposition  method  was  used  to  obtain  trends  of  the  six  cave  records.  After  the  base  trend  was  determined, 
81sO  anomalies  of  each  record  were  computed  by  subtracting  the  base  trend.  Mean  81sO  anomaly  values  of  the 
detrended  time  series  for  each  cave  record  were  calculated  for  four  periods:  (i)  Medieval  Warm  Period  (MWP;  AD  900- 
1 250):  (ii)  Little  Ice  Age  phase  1  (LIA-1 ;  AD  1 250-1  550);  (iii)  Little  Ice  Age  phase  2  (LIA-2;  AD  1  550-1 850);  and  (iv) 
modern  period  (MD;  AD  1  850-2000).  From  these  anomalies,  the  temporal  and  spatial  variability  of  wet/dry  conditions 
has  been  identified.  Positive  values  of  the  mean  81sO  anomalies  indicating  drier  conditions  appeared  in  the  lower 
Yangtze  River  Drainage  Area  and  Southeast  Coast  Area  during  MD-1,  LIA-1  and  MWP,  whereas  negative  values 
existed  in  north,  south  and  Yangtze  areas  of  eastern  China  during  LIA-2.  The  results  agree  with  the  dryness/wetness 
index  reconstructed  by  Chinese  historic  records  in  general.  These  results  illustrate  that  wet  and  dry  conditions  in 
different  regions  of  eastern  China  could  be  opposite  under  monsoon  influence,  so  that  no  single  speleothem 
81sO  record  could  represent  the  monsoonal  climate  in  this  vast  region.  Climatic  patterns  in  the  monsoonal  region 
can  be  either  a  combination  of  warm/wet  and  cold/dry  or  a  combination  of  cold/wet  and  warm/dry  on  annual-to- 
centennial  scales.  A  128-year  periodic  cycle  exists  in  all  six  cave  records,  whereas  64-year  (possibly  a  harmonic  of 
128-year  periodicity)  and  42-year  periodicities  appear  in  Shihua,  Heshang  and  Dongge  records.  These  cycles  may 
reflect  the  influence  of  solar  activity  on  the  East  Asian  Summer  Monsoon.  Copyright  ©  201 2  John  Wiley  &  Sons,  Ltd. 

KEYWORDS:  speleothem  8180  records;  East  Asian  Summer  Monsoon;  empirical  mode  decomposition  method;  Medieval  Warm 
Period;  Little  Ice  Age. 


Introduction 

The  East  Asian  Monsoon  (EAM)  is  the  annual  wind  direction 
change  caused  by  a  pressure  gradient  due  to  the  thermal 
contrast  between  ocean  and  land  surfaces.  These  pressure 
gradients  draw  moist  air  masses  landward  during  the  summer 
monsoon,  supplying  much  precipitation  to  eastern  China. 
During  the  winter,  flow  patterns  are  reversed,  and  cool,  dry  air 
dominates  over  eastern  China.  Long-term  fluctuations  in  East 
Asian  Summer  Monsoon  (EASM)  intensity  are  linked  to  glacial/ 
interglacial  cycles  (change  of  ice  volume)  and  orbital  forcing. 
Precisely  dated  speleothem  records  have  provided  a  history  of 
EASM  with  lighter  81sO  denoting  stronger  EASM  with  heavier 
precipitation  (Hu  etal.,  2008;  Ku  and  Li,  1998;  Li  etal.,  1998a, 
2011;  Paulsen  et  al.,  2003;  Wan  et  a!.,  2011a;  Wang  et  a!., 
2001;  Yuan  et  al.,  2004).  On  the  basis  of  extremely  sparse 
spatial  coverage  of  proxy  data  such  as  lacustrine,  marine 
sediments  and  speleothem  sequences,  it  is  known  that  the  EAM 
system  was  stronger  (weaker)  in  boreal  winters  but  weaker 
(stronger)  in  boreal  summers  during  glacial  (interglacial) 
periods  (e.g.  Ding  et  al .,  1995;  Wang  et  al.,  2008,  2010). 
For  example,  using  5-year-resolution  absolute-dated  spe¬ 
leothem  81sO  data  at  one  location  (Dongge  Cave)  in  southern 
China,  Wang  etal.  (2005)  found  pronounced  changes  in  EASM 
intensity  for  the  past  9000  years,  which  broadly  follows  the 
summer  insolation,  and  is  also  punctuated  by  eight  weak 
monsoon  events  lasting  around  one  to  five  centuries.  Up  until 
now,  almost  all  the  studies  on  speleothem  evidence  for  the 
EASM  change  are  concentrated  on  the  temporal  variation. 
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Zhang  et  al.  (2008)  used  81sO  in  stalagmite  WX42B  from 
Wanxiang  Cave  located  in  Gansu  Province  of  China  as  a  proxy 
for  EASM  strength  and  i  nterpreted  the  effect  of  changes  i  n  EASM 
strength  on  the  rise  and  fall  of  Chinese  dynasties.  However, 
monsoonal  rainfall  in  eastern  China  under  the  influence  of  the 
EAM  is  characterized  by  large  variability  in  space  and 
time  (Ding  and  Ren,  2008;  Guo  et  al.,  2003;  Shi  et  al., 
1 996).  The  influence  of  temporal  change  of  EASM  strength  has 
spatial  variation  (Zhang  et  al.,  2010,  and  references  therein). 
Figure  1  shows  the  correlations  of  June-August  rainfall  patterns 
with  summer  monsoon  strength  based  on  the  data  of  1 60 
meteorological  stations  in  China  (Guo  et  al.,  2003).  Summer 
monsoon  indices  (SMI)  were  reconstructed  using  sea-level 
pressure  difference  (<— 5  hPa)  between  land  (1 10°  E)  and  sea 
(160°  E)  over  latitudes  of  10-50°  N  in  June-August.  The 
influence  is  expressed  by  rainfall  during  the  10-year  period 
(1988-1997)  of  weakest  EASM  minus  rainfall  at  the  same 
station  during  the  10-year  period  (1955-1964)  of  strongest 
EASM.  It  is  seen  that  weaker  EASM  gave  rise  to  less  rainfall  in 
northern  China,  but  more  rainfall  in  the  middle-lower  reaches 
of  the  Yangtze  River  and  in  southern  coastal  areas.  Using 
dryness/wetness  index  (DWI)  reconstructed  from  historic 
records,  Zhang  et  al.  (2010)  also  obtained  spatial  variations 
of  the  monsoonal  rain  in  four  different  regions.  Therefore, 
modern  meteorological  observations  and  historic  records  tel  I  us 
that  a  strong  summer  monsoon  can  cause  wet  conditions  in  a 
region  but  dry  conditions  in  another  region  in  eastern  China. 
No  single  speleothem  81sO  record  can  represent  the  variation 
in  monsoonal  precipitation  of  the  whole  of  eastern  China. 

Recent  debates  about  the  physical  meaning  of  speleothem 
81sO  record  in  the  paleoclimate  community  have  focused  on 
how  speleothem  81sO  reflects  summer  monsoon  intensity  and 
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Figure  1.  Influence  of  EASM  on  summer  rainfall  in  China  during  AD 
1951-2000  (modified  from  Guo  ef  al.,  2003).  The  influence  is 
expressed  by  rainfall  during  the  10-year  period  (1 988-1 997)  of  weakest 
EASM  minus  rainfall  at  the  same  station  during  the  10-year  period 
(1955-1964)  of  strongest  EASM.  The  thick  contours  with  number  0 
indicate  no  difference.  The  thin  contour  lines  with  numbers  show 
rainfall  differences  (mm).  Crossed  areas  give  positive  values  and  dotted 
areas  negative  values.  Triangles  denote  locations  of  the  six  caves 
discussed  in  the  paper. 

rainfall.  As  changes  in  temperature,  rainfall,  upper  stream 
rainout  effect,  moisture  source  pathway  and  81sO  of  moisture 
source  all  can  affect  speleothem  81sO  record  (Wan  ef  al., 
2011a,  2011b,  and  references  therein),  a  different  frequency 
of  speleothem  81sO  record  may  have  different  main  forcing 
factors.  For  instance,  a  low  frequency  (>  1 03  years)  of 
speleothem  8180  record  variation  may  reflect  regional  or 
global  forcing  such  as  changes  in  global  ice  volume,  solar 
insolation  and/or  ocean  circulation,  whereas  high-frequency 
variations  in  the  speleothem  81sO  record  may  reflect  local 
rainfall  changes.  Thus  finding  an  appropriate  analytical  method 
for  understanding  the  different  frequency  bands,  anomaly  after 
detrending  and  cyclicity  of  a  speleothem  81sO  record  will  help 
us  to  explain  the  physical  meaning  of  the  speleothem 
81sO  record  in  general. 

To  what  extent  do  the  Chinese  speleothem  81sO  records 
provide  climatic  condition  in  terms  of  precipitation?  How  does 
one  extract  the  variation  of  monsoonal  rainfall  from  a 
speleothem  81sO?  Is  there  any  evident  spatial  variability  of 
EASM  in  short  duration  since  the  late  Holocene  represented  by 
Chinese  speleothem  records?  To  answer  these  questions,  high- 
resolution  81sO  data  (yearly)  since  AD  900  constructed  from 
speleothems  at  six  caves  (Dongge,  Furong,  Heshang,  Buddha, 
Shihua  and  Wanxiang)  in  eastern  China  are  used  to  identify  the 
spatial  inhomogeneity  of  EASM  fluctuation.  In  this  study,  we 
will  use  empirical  mode  decomposition  (EMD)  (Huang  et  al., 
1998)  to  detect  the  trends  of  the  six  high-resolution 
81sO  records.  The  anomalies  of  the  81sO  time  series  (relative 
to  the  trends)  are  determined  at  each  cave  and  the  dipole 
pattern  is  obtained  after  calculating  the  mean  anomalies  at  four 
climate  periods:  Medieval  Warm  Period  (MWP),  Little  Ice  Age 
phase  1  (LIA-1),  Little  Ice  Age  phase  2  (LIA-2)  and  modern 
period  (MD).  We  then  carry  out  spectral  analysis  on  the 
anomalous  81sO  records  to  detect  any  cycles  and  discuss  their 
significance.  Since  the  six  records  have  good  chronology 
control  and  the  age  uncertainties  of  the  records  are  much 
smaller  than  the  study  durations,  the  age  problem  is  not  an  issue 
here.  In  order  to  evaluate  our  results,  the  DWIs  reconstructed 
from  historic  documents  in  eastern  China  have  been  analyzed. 
Our  study  initiates  a  new  approach  to  using  high-resolution 


paleoclimatic  records  for  understanding  the  mechanisms  and 
forcing  factors  of  monsoonal  climate  changes. 

Speleothem  <S1sO  records  and  DWI 

Calcite  cave  deposits  can  effectively  record  climatic  variability 
occurring  at  different  timescales  over  the  Late  Pleistocene  (Bar- 
Matthews  et  al.,  2003;  Fairchild  et  al.,  2006;  McDermott, 
2004).  They  are  among  the  most  widely  distributed  in 
continental  environments,  and  are  amenable  to  precise  dating 
by  lamination  counting  and  uranium  series  (230Th/234U  and 
21  °Pb)  dating  (Baskaran  and  lliffe,  1993;  Edwards  ef  al.,  1986/ 
87;  Richards  and  Dorale,  2003;  Ludwig  et  al.,  1992;  Shen 
et  al.,  2002;  Shopov  et  al.,  1994;  Tan  et  al.,  2002).  The 
8180  fractionation  between  the  calcite  stalagmite  and  cave  drip 
water  is  a  function  of  cave  temperature  (7).  In  eastern  China, 
precipitation  mostly  (>75%)  occurs  in  the  summer  and  the 
vapor  source  is  mainly  the  western  North  Pacific  and  Indian 
oceans.  This  seasonality  and  moisture  source  of  the  monsoon 
system  have  remained  little  changed  during  the  late  Holocene. 
Therefore,  although  multiple  factors  such  as  temperature  and 
changes  in  moisture  source  and  its  8180  can  affect  the  81sO  of  a 
stalagmite,  rainfall  amount  is  the  dominant  factor  influencing 
stalagmite  81sO  on  annual  to  decadal  scales  in  eastern  China 
(Li  et  al.,  1 998a). 

High-resolution  (annual-to-decadal),  precisely  dated  spe¬ 
leothem  8180  records  of  the  six  caves  in  the  monsoonal  region 
of  eastern  China  were  selected  to  study  the  temporal-spatial 
variability  of  EASM  since  the  MWP  (Fig.  2).  Table  1  shows  the 
locations  of  these  caves  and  nearby  cities  where  meteorolo¬ 
gical  data  and  historic  climate  records  can  be  obtained.  All  six 
81sO  records  have  been  published:  the  Dongge  Cave  record  by 
Wang  et  al.  (2005),  the  Furong  Cave  record  by  Li  et  al.  (201 1 ), 
the  Heshang  Cave  record  by  Hu  etal.  (2008),  the  Buddha  Cave 
record  by  Paulsen  et  al.  (2003),  the  Wanxiang  Cave  record  by 
Zhang  et  al.  (2008)  and  the  Shihua  Cave  record  by  Wan  et  al. 
(201 1  b).  The  stalagmites  were  dated  by  various  methods  such 
as  230Th/U,  210Pb  and  lamination  counting  to  a  time  resolution 
as  fine  as  year-to-subdecade  for  the  last  millennium.  These 
records  were  generally  compared  with  local  rainfall  records, 
and  some  of  them  were  used  to  interpret  historic  events  (e.g.  Li 
ef  al.,  1998b;  Zhang  et  al.,  2008)  and  matched  with  solar 
variability  (e.g.  Wang  etal.,  2005).  Figure  3  (light  curves)  shows 
the  original  speleothem  81sO  records.  Yearly  speleothem  8180 
time  series  of  the  six  records  were  obtained  through  a  linear 
interpolation  between  points  that  are  not  annual  resolution.  In 
this  study,  time  durations  of  the  MWP,  early  and  late  LIA  and 
modern  since  the  Industrial  Revolution  were  chosen  (Crowley, 
2000;  Hughes  and  Diaz,  1 994;  Jansen  ef  al.,  2007).  To  work  on 
multi-century  durations,  the  chronologies  of  the  records  are 
precise  enough  within  the  uncertainties  (<±15  years). 

In  order  to  cross-check  our  results  obtained  from  the 
speleothem  81sO  records,  DWI  in  16  locations  as  well  as 
the  mean  DWI  in  the  three  divisions  of  Huabei,  Huanan  and 
Low-Middle  Yangtze  River  Drainage  Basin  were  selected 
(Table  1).  Chinese  historical  documentations  are  generally 
accurate  in  chronology  and  unambiguous  in  their  description  of 
dry/wet  and  warm/cold  conditions  (Zhang,  1988;  Zheng  et  al., 
1977,  2006).  These  records  are  based  on  abundant  historic 
documents  in  China  and  have  been  calibrated  with  instru¬ 
mental  precipitation  records.  In  this  study,  only  the  records 
back  in  AD  1471  or  beyond  were  selected  (Chinese  Academy 
of  Meteorological  Sciences,  1 981 ;  Zhang  ef  al.,  2003).  For  the 
DWI  records,  climatic  conditions  are  classified  into  five 
categories:  2  and  —2  stand  for  very  wet  and  very  dry, 
respectively;  1  and  —1  denote  for  wet  and  dry,  respectively;  and 
0  represents  normal. 
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Figure  2.  Location  map  showing  the  six 
caves  discussed  in  the  paper  that  have  high- 
resolution  8180  records  with  good  chronology 
control  and  1  6  stations  which  have  DWI  since 
AD  1471,  except  Beijing  record,  which  goes 
back  to  AD  960. 


Trends  of  speleothem  51sO  data 

To  identify  EASM  fluctuations,  the  trends  of  speleothem 
81sO  should  be  first  determined  for  the  six  cave  records.  Let 
x(t)  represent  the  time  series  of  the  speleothem  81sO  yearly  data 


as  shown  in  Fig.  3.  From  AD  900  to  AD  2000,  x(t)  fluctuates  on 
various  timescales  (inter-annual,  decadal  and  centennial)  with 
many  local  maxima  and  local  minima.  Huang  el  al.  (1998) 
developed  empirical  mode  decomposition  (EMD)  to  objec¬ 
tively  obtain  the  trend.  EMD  is  a  non-parametric  data-driven 


Table  1.  Information  of  six  caves  (*)  and  stations  for  historic  DWI. 


MD 

LIA-2 

LIA-1 

MWP 

Site  name 

Location 

2000-1 850 

1849-1550 

1549-1250 

1 250-960 

Huabei  Division 

-0.028 

0.076 

Beijing  (BJ) 

39.924° 

N,  116.381°  E 

-0.126 

0.076 

0.129 

0.004 

Tianjin  (TJ) 

39.131° 

N,  117.203°  E 

0.000 

0.130 

Shihua  Cave* 

39°  48' 

N,  115°  54'  E 

-0.126 

-0.014 

-0.006 

0.032 

Baoding  (BD) 

38.857° 

N,  115.500°  E 

0.046 

-0.050 

Low-Middle  Yangtze  River  Division 

0.099 

0.139 

Shanghai  (SH) 

31.238° 

N,  121.469°  E 

0.020 

0.153 

Suzhou  (SZ) 

31.316° 

N,  120.619°  E 

0.020 

0.143 

Yichang  (YC) 

30.704° 

N,  1  1  1.285°  E 

0.099 

0.147 

Heshang  Cave* 

30°  27' 

N,  110°  25'  E 

0.166 

-0.039 

0.001 

0.078 

Huanan  Division 

-0.025 

0.093 

Guilin  (GL) 

25.282° 

N,  110.287°  E 

0.053 

0.070 

Guangzhou  (GZ) 

Other  eastern  China  sites 

23.119° 

N,  113.261°  E 

-0.073 

0.093 

Datong  (DT) 

40.097° 

N,  113.296°  E 

0.066 

-0.101 

Taiyuan  (TY) 

37.871° 

N,  1  12.569°  E 

-0.026 

-0.013 

Linfen  (LF) 

36.083° 

N,  111.514°  E 

-0.086 

-0.110 

Changzhi  (CZ) 

36.182° 

N,  113.106°  E 

-0.205 

-0.110 

Yulin  (YL) 

38.297° 

N,  109.757°  E 

-0.159 

-0.193 

YanAn  (YA) 

36.594° 

N,  109.471°  E 

-0.166 

-0.167 

Xian  (XA) 

34.262° 

N,  108.949°  E 

0.060 

-0.067 

Buddha  Cave* 

33°  40' 

N,  109°  05'  E 

0.343 

-0.191 

0.013 

0.059 

Ankang  (AK) 

32.690° 

N,  109.026°  E 

0.278 

0.217 

Hanzhong  (HZ) 

33.078° 

N,  107.034°  E 

0.179 

0.233 

Wanxiang  Cave* 

33°  19' 

N,  105°  00'  E 

-0.053 

0.010 

-0.018 

-0.021 

Guiyang  (GY) 

26.577° 

N,  106.711°  E 

-0.106 

-0.013 

Furong  Cave* 

29°  13' 

N,  107°  54'E 

-0.208 

0.083 

0.032 

-0.004 

Dongge  Cave* 

25°  1  7'N,  108°  50'E 

-0.043 

-0.013 

0.013 

-0.042 

All  locations  of  the  caves  and  stations  are  shown  in  Fig.  2.  Three  divisions  including  Huabei  (northern  China),  Huanan  (southern  China)  and  Low- 
Middle  Yangtze  River  Drainage  Basin  are  climatic  regions  based  on  the  meteorological  observation  over  the  past  60  years  in  China.  The  other  sites  also 
belong  to  monsoonal  area  of  eastern  China,  but  outside  of  the  above  three  divisions.  The  calculated  mean  climatic  conditions  for  the  four  periods  are 
shown  by  the  numerical  numbers.  For  DWI,  positive  numbers  indicate  wet  climate  and  negative  numbers  are  dry  climate.  For  the  cave  8180  anomaly, 
positive  values  reflect  dry  climate,  whereas  negative  numbers  show  wet  climate. 
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Figure  3.  Time  series  of  high-resolution 
(yearly)  speleothem  8180  data  from  south  to 
north  at  the  (a)  Dongge,  (b)  Furong,  (c) 
FHeshang,  (d)  Buddha,  (e)  Wanxiang  and  (f) 
Shihua  caves.  FHere,  the  dashed  curves 
represent  the  trends  obtained  using  the  EMD 
method.  It  is  noted  that  all  the  trends  show  a 
similar  pattern  (except  the  recent  increase  of 
speleothem  8180  since  1800  except  at  the 
Wanxiang  cave). 


(a)  Dongge 


Year 


analysis  tool  that  decomposes  nonlinear  non-stationary  signals 
into  intrinsic  mode  functions  (IMFs).  An  IMF  is  a  function  that 
must  satisfy  two  conditions  according  to  the  algorithm 
originally  developed:  (a)  the  difference  between  the  number 
of  local  extrema  and  number  of  zero  crossings  must  be  zero  or 
one;  (b)  the  running  mean  value  of  the  envelope  defined  by  the 
local  maxima  and  the  envelope  defined  by  the  local  minima  is 
zero.  The  algorithm  to  decompose  a  signal  into  IMFs  is  then  the 
following  (FHuang  et  a!.,  1998):  First,  the  local  minima  and 
maxima  of  the  signal  x{t)  are  identified.  Second,  the  local 
maxima  are  connected  together  by  a  cubic  spline  interpolation 
(other  interpolations  are  also  possible),  forming  an  upper 
envelope  emax(f).  The  same  is  done  for  local  minima,  providing 
a  lower  envelope  emin(f).  Third,  the  mean  of  the  two  envelopes 
is  calculated: 

(0  =  [^max(0  +  (1 ) 

Such  a  procedure  is  shown  in  Fig.  4.  Fourth,  the  mean  is 
subtracted  from  the  signal,  providing  the  local  detail: 

/?!  (0  =  x(f)  -  (t)  (2) 


&nra»opK  and  Ihn  Mann ;  tfala 


Figure  4.  Illustration  of  the  EMD  process,  with  the  blue  curve 
denoting  x(f),  the  two  green  curves  representing  envelopes  of  local 
maxima  emax(f)  and  minima  emin(f)  and  the  red  curve  referring  to 
m,(r)=  [emax(f)  +  emin(t)]/2. 
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which  is  then  considered  to  check  whether  it  satisfies  the  above 
two  conditions  to  be  an  IMF.  If  yes,  it  is  considered  as  the  first 
IMF  and  denoted 


X,  (f)  =  /i,  (f) 


(3) 


It  is  subtracted  from  the  original  signal  and  the  first  residual: 

n  (t)  =  x(f)  -  X\  (f)  (4) 

which  is  taken  as  the  new  series  in  step  1 .  If  A?n  (f)  is  not  an  IMF,  a 
procedure  called  a  'sifting  process'  is  applied  as  many  times  as 


(d)  Buddha  Cave 


Year 


Year 


^  ^  A  ^ 

Year 


Figure  5.  Functions  c-|(f)  to  c7(f)  and  the  trends  of  the  speleothem  81sO  time  series  at  the  (a)  Dongge,  (b)  Furong,  (c)  Heshang,  (d)  Buddha,  (e) 
Wanxiang  and  (f)  Shihua  caves. 
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necessary  to  obtain  an  IMF.  In  the  sifting  process,  h^(t)  is 
considered  as  the  new  data,  and  the  same  procedure  applies. 
The  IMFs  are  orthogonal,  or  almost  orthogonal  functions 
(mutually  uncorrelated).  This  method  does  not  require 
stationarity  of  the  data  and  is  especially  suitable  for  non¬ 
stationary  and  nonlinear  time  series  analysis. 

By  construction,  the  number  of  extrema  decreases  when 
going  from  one  residual  to  the  next;  the  above  algorithm  ends 
when  the  residual  has  only  one  extrema,  or  is  constant,  and  in 
this  case  no  more  IMF  can  be  extracted;  the  complete 
decomposition  is  then  achieved  in  a  finite  number  of  steps. 
The  signal  x(t)  is  finally  written  as  the  sum  of  mode  time  series 
Xj(t)  and  the  residual  rm(t): 

m 

x{t)  =  J2xi(t)  +  rm(t)  (5) 

1=1 

where  x-i(f)  has  the  highest  temporal  variability  and  xm(t)  has  the 
lowest  temporal  variability.  The  functions 

k 

Ck(t)  x(f)  -  y>(f),  k  =  1,2,..., m  (6) 

/= i 

show  the  filtration  of  high-frequency  variability  from  the  signal 
x(t)  with  c-i(t)  filtering  out  of  x-| (t),  c2(f)  filtering  out  of 
x1(f)  +  x2(f),  ...,  and  cm(t)  is  rm(t). 

A  residual  rm(t)  can  be  treated  as  a  trend  if  its  low-frequency 
variability  occurs  in  all  the  functions  q(f),  c2(f),  . . .,  cm{t).  In 
other  words,  the  trend  should  not  be  filtered  out  by  any 
methods.  Figure  5  shows  such  a  filtration  process  using  the 
EMD  method  on  the  speleothem  81sO  time  series  at  the  six 
caves.  It  is  noted  that  the  low-frequency  variability  represented 
by  r7{t)  (identified  by  the  trend)  occurs  in  c-|(f),  c2(f),  . . .,  c/t). 

To  identify  the  temporal  variation  in  long  timescales,  four 
major  periods  are  defined:  MWP,  LIA-1 ,  LIA-2  and  MD.  Trends 
of  the  speleothem  81sO  time  series  are  similar  in  these  caves 
(except  for  the  Wanxiang  Cave)  (dashed  curves  in  Fig.  3):  it 
decreases  slightly  with  time  during  MWP  (AD  900-1 250),  with 
minimum  values  occurring  around  AD  1 1  50  (strengthening  of 
the  EASM);  increases  with  time  during  LIA-1  (AD  1 250-1  550), 
with  maximum  values  appearing  around  AD  1 500  (weakening 
of  the  EASM);  decreases  with  time  during  LIA-2  (AD  1550- 
1850),  with  minimum  values  occurring  around  AD  1800 
(strengthening  of  the  EASM);  and  increases  with  time  during 
MD  (AD  1850-2000)  (weakening  of  the  EASM).  FHowever,  the 
trends  of  the  speleothem  81sO  time  series  for  the  Wanxiang 
Cave  show  an  evident  decrease  during  MD  (Fig.  3). 

Speleothem  51sO  anomaly 

The  detrended  data, 

x(f)  =  x(f)  -  rm(t)  (7) 

is  called  the  speleothem  81sO  anomaly,  which  changes  on 
various  timescales  at  the  six  caves  (Fig.  6).  With  these  time 
series  of  speleothem  81sO  anomaly,  spectral  analysis  has  been 
conducted.  The  power  spectra  with  a  confidence  level  of  95% 
show  an  evident  period  of  1 28  years  for  all  six  caves,  and  other 
periods  such  as  42-year  and  64-year  cycles  in  different  records 
(Fig.  7). 

Many  studies  of  natural  archives  have  shown  the  wide 
periodic  dominance  in  the  40-  to  128-year  band,  e.g.  in  tree 
rings  (Gray  et  al.,  2004;  Ogurtsov  et  al.,  2002;  Ware  and 
Thomson,  2000),  estuarine  fossil  pigments  (Hubeny  et  al., 
2006),  marine  sediments  and  fish  scales  (Berger  et  al.,  2004).  In 
tree  ring  records,  the  42-year  and  64-year  periodicities  were 
found  to  reflect  the  Atlantic  Multidecadal  Oscillation  (AMO), 
which  plays  a  role  in  the  North  Atlantic  climate  by  producing 


Figure  6.  Detrended  speleothem  8180  time  series  at  the  six  caves. 


anomalous  geopotential  heights  over  the  fall  and  summer 
seasons  (Enfield  and  Mestas-Nunez,  1999;  Gray  et  al.,  2004; 
FHubeny  et  al.,  2006).  Similar  to  the  AMO,  multidecadal 
variability  of  climatic  proxies  is  found  in  western  North 
America  (Gray  et  al.,  2003)  and  the  Northwest  Pacific  Basin 
(Delworth  and  Mann,  2000;  Ware  and  Thomson,  2000).  It  is 
generally  accepted  that  the  64-  to  140-year  band  may  reflect 
the  respective  changes  in  the  length  of  the  solar  Gleissberg 
cycle,  which  is  often  centered  at  88-98  years  (Gleissberg, 
1944;  Ogurtsov  et  al.,  2002).  Sunspot  numbers  (e.g.  Wolf 
numbers)  have  wide  periodic  dominance  in  the  64-  to  1 28-year 
band  by  wavelet  filtered  analysis  (Kane,  1 999;  Ogurtsov  et  al., 
2002).  Despite  some  debate  about  whether  the  small  variation 
in  solar  activity  could  significantly  affect  the  low  atmospheric 
climates,  many  studies  show  that  the  multidecadal  variability  of 
regional  climates  may  be  attributed  to  the  solar-oceanic- 
atmospheric  connections.  For  instance,  multidecadal  oceanic 
fluctuations  associated  with  the  AMO  are  teleconnected  with 
Pacific  Basin  modes  via  the  AO  and  a  hypothesized 
atmospheric  bridge  at  high  latitudes  (FHonda  et  al.,  2001).  In 
combination  with  tropical  Pacific-lndian  Ocean  forcing  and 
internal  processes,  the  AMO  is  influenced  by  the  North  Atlantic 
Oscillation  (NAO)  and  Arctic  Oscillation  (AO),  which  are 
partially  driven  by  North  Atlantic  sea-surface  temperatures 
(SSTs)  (FHoerling  et  al.,  2001;  Sutton  and  FHodson,  2003).  It  is 
acceptable  that  changes  in  solar  activity  such  as  sunspot 
numbers  can  cause  variations  of  SST,  which  further  affect 
atmospheric  circulations  by  altering  evaporation,  precipitation 
and  ocean-atmosphere  heat  exchanges.  Delworth  and  Mann 
(2000)  showed  that  the  AMO  is  related  to  changes  in 
thermohaline  circulation  and  meridional  heat  flux.  FHence 
we  may  hypothesize  that  changes  in  solar  activity  may  affect 
SST  and  atmospheric  circulations  on  ocean  basins  such  as 
NAO,  Pacific  Decadal  Oscillation  (PDO)  and  AO.  The  surface 
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Spectrum  Analysis 
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Figure  7.  Power  spectra  of  the  detrended  speleothem  8180  time  series  at  the  six  caves.  Here  the  horizontal  axis  shows  the  period.  It  is  noted  that  the 
period  of  128  years  is  very  evident  in  all  the  caves. 


changes  lead  to  changes  in  thermohaline  circulation  and 
meridional  heat  flux  that  cause  multidecadal  climate  fluctu¬ 
ations. 

It  is  interesting  to  see  that  all  six  records  have  128-year 
periodicity  (Fig.  7).  This  lower  frequency  band  may  imply  that 
the  solar  variability  has  a  broad  impact  on  the  EASM,  with 
minor  modification/interruption  of  local  or  regional  complexes 
in  climatic  settings.  If  considering  the  64-year  periodicity  as  a 
harmonic  of  1 28-year  periodicity,  the  64-year  periodicity  could 
appear  in  all  six  records.  However,  this  64-year  cycle  is  not 
shown  in  the  Wanxiang  and  Furong  Cave  records.  Actually,  the 
shorter  periodicities  (<1 28  years)  are  not  seen  in  Wanxiang  and 
Furong  Cave  records,  which  are  the  two  further  inland  cave 


sites  (Fig.  7).  If  this  is  the  actual  situation,  it  may  imply  that 
the  shorter  periodicities  (64  and  42  years)  are  probably  more 
related  to  oceanic  changes  which  affect  the  climate  of  the  area 
closer  to  the  ocean  (e.g.  Shihua  and  Dongge  Cave  sites).  In 
view  of  the  physical  meanings  of  speleothem  81sO,  variation 
of  the  81sO  on  longer  timescales  may  not  reflect  rainfall 
amount  change  but  moisture  source,  etc.,  under  the  influence 
of  the  summer  monsoon.  The  128-year  cycle  thus  may 
indicate  the  speleothem  81sO  change  caused  by  changes  in 
moisture  source  (including  its  81sO  and  pathway,  such  as 
Pacific  vs.  Indian  Ocean).  On  this  timescale,  the  speleothem 
81sO  may  not  present  rainfall  amount  change.  In  this  case,  one 
should  think  about  the  climatic  implication  of  the  similarity  of 
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the  Holocene  trend  in  most  speleothem  81sO  records.  The  high 
frequent  fluctuations  of  the  speleothem  81sO  may  reflect  mainly 
local  rainfall  changes.  Figure  7  shows  that  the  64-year  and 
42.7-year  cycles  are  clear  monsoonal  rainfall  cycles  in  Shihua, 
Dongge  and  Heshang  Cave  sites,  which  are  closer  to  the  ocean. 
The  25.6-year  cycle  in  the  Shihua  Cave  record  (northernmost) 
and  32-year  cycle  in  the  Dongge  Cave  record  (southernmost) 
show  more  regional  differences  in  rainfall  cycles  even  though  the 
two  sites  are  close  to  the  coastal  area. 

To  identify  the  temporal  variation  in  the  four  major 
periods  (MWP,  LIA-1,  LIA-2  and  MD),  the  speleothem 
81sO  anomaly  x(t)  is  averaged  within  the  four  periods 
(*mwp,*lia-i,*lia-2,*md)  for  each  cave  (Table  1).  Figure  8 
shows  the  spatial  distribution  of  xMwp5 Xlia-i, Xlia-2, xMd/ 
respectively. 

The  speleothem  81sO  anomaly  (or  the  monsoonal  rainfall 
anomaly)  reveals  an  east-west  dipole  pattern  with  alternating 
(+,  — )  signs  as  the  climate  period  shifts.  During  the  MWP 
and  LIA-1,  the  eastern  part  of  eastern  China  has  a  positive 
anomaly  (xmwp  >  0,  xLia-i  >  0),  referring  to  drier  climatic 
condition;  and  the  western  part  of  eastern  China  has  a  negative 
anomaly  (xMwp  <  0,  Xlia-i  <  0),  referring  to  wetter  conditions. 
During  the  LIA-2,  the  pattern  flips,  in  that  the  eastern  part  of 
eastern  China  has  negative  anomaly  (xLIA-i  <  0),  referring  to 
increase  of  summer  rainfall;  and  the  western  part  of  eastern 
China  has  a  positive  anomaly  (xLia-i  >  0),  referring  to  decrease 
of  summer  rainfall.  During  the  MD,  the  pattern  flips  back 
(similar  to  the  MWP  and  LIA-1),  in  that  the  eastern  part  of 


eastern  China  has  a  positive  anomaly  (xMd  >  0);  and  the 
western  part  of  eastern  China  has  a  negative  anomaly 
(*md  <  0). 

For  DWI  records,  the  mean  values  for  each  station  at  four 
periods  are  calculated  and  listed  in  Table  1 .  Only  Beijing  DWI 
is  1000  years  long,  so  that  the  rest  of  the  DWI  records  have  only 
mean  values  for  periods  of  LIA-2  and  MD.  For  LIA-2,  which  is 
the  second  half  of  the  Little  Ice  Age,  the  mean  DWI  values  in  all 
three  divisions  are  positive,  with  mean  a  wet  climate  in  this 
period  in  the  eastern  part  of  eastern  China.  This  agrees  with  the 
speleothem  81sO  records,  which  showed  a  negative  value  for 
wetter  climatic  conditions  during  LIA-2.  In  contrast,  the  mean 
DWI  values  of  the  most  western  sites  (DT,  TY,  LF,  CZ,  YL,  YA, 
XA,  GY)  were  negative  numbers,  showing  dry  climates  in  the 
same  period.  We  admit  that  the  cave  records  are  too  sparse, 
especially  for  the  northern  part,  so  that  the  curves  drawn  in 
Fig.  8  are  not  accurate  enough  to  compare  with  the  DWI 
records.  For  LIA-1,  the  Beijing  DWI  record  supports  the 
speleothem  reconstruction,  showing  a  wet  climate  in  that  area. 
However,  for  the  MWP  the  cave  record  indicates  a  dry 
condition  in  Beijing  area  but  the  Beijing  DWI  shows  a  small 
positive  value,  which  may  have  some  uncertainty.  Never¬ 
theless,  the  Beijing  mean  DWI  value  does  not  show  a  significant 
dry  climate  during  MWP.  The  comparison  of  the  two  kinds 
of  records  for  the  modern  period  (MD)  do  not  compare  well. 
The  mean  DWI  values  of  most  sites  in  eastern  China  are 
negative  numbers  except  for  the  Low-Middle  Yangtze  River 
Division,  reflecting  dry  climate  conditions.  Mean  values  of  the 


(a)  MWP  (AD1125-AD900)  (b)  LIA(phase-l)  (AD1550-AD1125) 


(c)  LIA(phase-ll)  (AD1850-AD1550)  Modern  (AD2000-AD1850) 


Longitude 


Longitude 


Figure  8.  Spatial  distributions  of  mean  81sO  anomalies  during  (a)  MWP,  (b)  LIA-1,  (c)  LIA-2  and  (d)  MD. 
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speleothem  81sO  anomalies  in  the  western  part  of  eastern 
China  are  negative,  referring  to  wetter  conditions.  Although  the 
global  warming  trend  and  industrial  pollution  can  affect  both 
the  relationship  of  speleothem  8180-temperature  and  the 
relationship  of  monsoonal  strength  with  rainfall  pattern,  the 
spatial  variation  of  the  EASM  strength  during  the  MD  shown  in 
Fig.  8  needs  more  evidence  for  clarification. 

Although  the  distribution  patterns  of  the  EASM  climates 
during  the  four  periods  are  roughly  reconstructed,  it  is  clear  that 
the  summer  monsoon  strength  cannot  produce  a  similar 
moisture  condition  over  eastern  China  on  annual-to-centennial 
timescales.  The  spatial  patterns  also  change  with  time.  No 
single  speleothem  record  can  interpret  the  paleoclimate  for  the 
whole  of  eastern  China.  In  addition,  the  monsoonal  climatic 
patterns  are  not  only  warm/wet  and  cold/dry,  but  also  can  be 
warm/dry  and  cold/wet,  depending  on  location  and  time.  In 
order  to  reconstruct  detailed  monsoonal  climates  and  to 
understand  their  forcing  factors,  more  high-resolution  spe¬ 
leothem  records  distributed  in  broad  regions  are  needed. 

Conclusions 

Temporal  and  spatial  variability  of  the  EASM  climate  (in  terms 
of  rainfall)  since  AD  900  has  been  identified  from  high- 
resolution  speleothem  8180  data  (yearly)  at  six  caves  (Dongge, 
Furong,  Heshang,  Wanxiang,  Buddha,  Shihua)  in  eastern 
China.  The  trends  of  the  speleothem  8180  time  series  at  these 
caves  were  obtained  using  the  EMD  method.  The  trends  of 
speleothem  81sO  imply  strengthening  of  the  EASM  during  the 
MWP  (AD  900-1 250)  and  LIA-2  (AD  1 550-1 850),  respectively; 
and  weakening  of  the  EASM  during  LIA-1  (AD  1 250-1  550)  and 
MD  (AD  1850-2000),  respectively.  The  detrended  high- 
resolution  speleothem  8180  data  (i.e.  anomaly)  show  local 
climatic  conditions,  with  a  negative  anomaly  reflecting  wet 
climates  and  vice  versa.  The  spectral  analysis  on  the 
speleothem  81sO  anomaly  time  series  shows  that  all  six 
records  have  a  128-year  periodic  cycle,  and  three  locations 
close  to  the  ocean  have  64-year  and  42-year  periodicities  as 
well.  These  cycles  are  probably  evidence  of  solar  influence  on 
the  EASM.  The  temporally  averaged  speleothem  81sO  anomaly 
for  the  four  climate  periods  during  the  last  millennium  reveals 
an  east-west  dipole  pattern,  with  alternating  (+,  — )  signs  as  the 
climate  period  shifts:  dry  climates  (wet  climates)  in  the  eastern 
(western)  part  of  eastern  China  during  the  MWP  and  LIA-1 ,  wet 
climates  (dry  climates)  in  the  eastern  (western)  part  of  eastern 
China  during  the  LIA-2,  and  dry  climates  (wet  climates)  in  the 
eastern  (western)  part  of  eastern  China  during  the  MD.  Therefore 
the  monsoonal  climatic  patterns  have  different  combinations  in 
terms  of  temperature  and  wetness  depending  on  location  and 
time.  No  single  speleothem  record  can  represent  the  climatic 
conditions  in  a  broad  region  of  eastern  China. 
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